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Abstract 

We calculate the next-to-leading corrections to the branching ratio of exclusive B — > 
A"*7 decay. The renormalization scale dependence is reduced compared to the leading 
logarithmic result but there remains a dependence on a cutoff parameter of the hadronic 
model. The calculated corrections increase the predicted branching ratio by about 10%, 
but it remains in agreement with the experimental value. 



1 Introduction 



The rare neutral flavor changing decays of B-meson are generally considered a good testing 
ground for new physics The inclusive decay b — > S7 has therefore received intensive 
experimental || and theoretical attention p||- 1|1 1|| . In particular, the next-to-leading order 
calculation for the standard model strongly improved the theoretical uncertainties, most 
notably the scale dependence of the perturbative calculation. 

The calculation of the decay rates for the corresponding exclusive decays such as 
B — > K*^( Jl^] is usually much more complex because of bound state effects (for a recent 



work, see for instance |I3| ). But they are easier to observe, and so it is also of interest 
to calculate their branching ratios and direct CP asymmetries to the highest possible 
precision. In this paper we consider the branching ratio for B — > K*^ in the next-to- 
leading logarithmic approximation. The decay B — > p7 can in principle be treated in the 
same way, but there are several additional contributions which require new methods. 



In ||14|| , the B — > /ry and B — > K*^ branching ratios and rate asymmetries were 
investigated in the leading logarithmic approximation within a simple meson model. Al- 
though this model does not give completely correct results (in particular the ratio of the 
semileptonic tt and p decays), we will also use it in this work, because of its simplicity. 

As usual, the basis of the calculation is the effective Hamiltonian which we record here 



for completeness. For B — > K*j decay it has the form [14 



Heff = {ZiVtdMOiQj,) - v u [dfaXOM - OM) + C 2 (fi)(0 2 M - 2 (ja))]} , 

(1) 

where v a = V t bV* s , a=u,c,t, i=l,2,...8. For the decay B — > p7 one must replace v a by v a , 
with v a = V t bV* d . The operators Oi can be found in many places |8|, [1^, |16| . 



2 The bound state model 

As mentioned, we will use a simple model for the mesons to take into account the bound 
state effects [0, [17| (for other approaches, see |1|,[|19|). The B and the K* are described 
by two constituents: B = (b,q), K* = (s,q). The sum of the four momenta of the 
constituents equals that of the bound state. The mass of the spectator antiquarks is 
taken to be zero and the b-quark mass becomes momentum dependent: 



m b = \Jmr B - 2p B p q , (2) 

where p B and p q are the four momenta of the B-meson and spectator quark, respectively. 
The momentum p q in the B-meson rest frame is restricted to 

\p q \ =E< m B /2. (3) 

The B meson is represented by a matrix ty B 



( / TTiS^ x I ^T^ZTZ-^b ( \p B Pg I /m B ) x -^=1 



2E,(2tt)3V ^{ml-pBVa) Kl qU VN' 
S B = ~(pb + m b )>y 5 pg, <f) B (p) = exp(-p 2 /(2p 2 F )), (4) 
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where C B is a normalization factor and pp is fixed in such a way that one obtains the 
correct value of B- meson decay constant f B . The 'hat' symbol stands for the usual 'slash'. 

The final state vector meson K* is represented by two constituents with parallel mo- 
menta p s = ypv, P q = (1 — y)pv- 111 most cases its mass and that of the constituent 
s quark can be neglected (see however the discussion of that point in section 4). It is 
described by the matrix 

■Vy = C v f 1 dy<f> v (y)X v x -^LlAr, C v = f v /{WN) 
Jo JN 



£ y = lim Jv{Pv + m v ), <f> v (y) = 6y(l - j/)(H ), (5) 

where 4>v(y) for V = K* , p are given by [[0| |20[| : 

<fo.(y) = 6y(l- 2 /)(l + 0.57(2y-l)-1.35((2 2 /-l) 2 -l/5) + 

+ 0.46(7(2?/ - l) 3 /3 - (22/ - 1))) (6) 
Mv) = 62/(1 -y)(l-0.85((2y-l) 2 - 1/5)). 
The B — > K*j transition matrix elements have a following general form: 

< Vj\Oi\B >= em B e^ [ie^p^eyPyF^Oi] + ((p 7 e v )p Vlx - (p 7 p y )ey M )F 5 [O i ]] . (7) 

For a massless final state meson we have F\ — F$. 

In our model, the hadronic matrix elements are written as 

<v imB > = CbCv /^^/=^= 

J [2'nY2E q y 2(m B - p B p q ) 

x <t> B (^)My)^Tr[Z v M sb Z B M qq , 

where the trace is in Dirac and color space and the matrices M sb and M qq i are related to 
the quark-level matrix elements by 

< sqj\Oi\bq >= u s M sb u b v q M qq >v q > . (9) 

The decay width in NLL approximation can be written as 

T (B - Vi) = \v t \ 2 G 2 F a em m B £ | £CiF,[0 4 ]| a (l - . (10) 

In leading approximation only Oj contributes and the last parenthesis is simply one. 
In order to eliminate some of the uncertainties of the model, we write 

BR(B - V-y) = T{B ^ Vl) BR(B -> Xlvi), (11) 

J- si 

where T s i in NLL approximation is given by 

g (r) = 1 - 8r 2 + 8r 6 - r 8 - 24r log(r) (12) 

3 2)(l_ r )2 + 3. 

4 ; 2' 



si 



f(r) = (7r 2 - — )(1 - r) 2 + -, r = m c /m b (p q ) 



and we use the experimental value for the semileptonic decay, BR(B — ► X/i/k)=0.103 |21 
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3 The "inclusive amplitudes" 

Consider first the decay graphs where the spectator is not touched; they are very similar 
to the inclusive decays. As mentioned, in leading logarithmic approximation only F[0 7 ] 
contributes to < Vj\T-C\B > (Fig. 1). Evaluating the trace, one gets JT3J]: 

[ ° 7] = J dy^f^My)^ ^ g j ■ (13) 

Here, and everywhere we will set F\ = F 5 = F for all form factors. The deviation due to 
small final state masses can be neglected. 

For the NLL amplitudes, we can use the results for NLL inclusive 7 , 2 , 2 ampli- 
tudes given in || 

< S1 \0 7 \b >=< S1 \0 7 \b > tree (l + ( 7 }° 7 )e// Iog(WA*) + rrfj 

< s^\0 2 \b >=< s-/\0 7 \b >tree ^J 12 ^ (l 2 \og{m b / 'n) + r 2 ) 

< S1 \0 8 \b >=< S1 \0 7 \b > tree -^-^(l 8 log{m b /fi)+r 8 ) (14) 

47T 

< s-i\0 7 \b >tree= m b ^—u{p')equ(p), 

where the expressions for r 2 , r 7 , r 8 , l 2 , , y 77 e ^ , h are given in ||. 

The problem of Bremsstrahlung corrections needs discussion. In the inclusive case the 
infrared divergences from the NLL calculation (loops) of < S7IO7I& > (Fig. 2) cancels at 
the level of the decay width when one adds Bremsstrahlung corrections associated with 
the operator 7 (that is when we sum over final states with additional soft gluons). In the 
exclusive case the situation is much more complicated. Normally we consider the mesons 
to contain only two constituents (quark and antiquark) in a color singlet state. Therefore, 
it is not possible to simply add to the Hilbert space of hadrons states with extra gluons; 
furthermore, soft gluons are not physical. Instead, these gluons must be incorporated 
into the physical hadrons, making it necessary to consider three particle states (quark, 
antiquark and gluon) etc for the K* meson and other physical particles. For a discussion 
of this point see p2| . 

This well known problem cannot be really solved with the present methods and one 
must resort to approximate methods (Additional IR problems connected to the exchange 
of gluons with the spectator quark will be discussed shortly). 

We have therefore used two heuristic methods to obtain finite results. In Method 1 
we use the fact that we could write the finite 'total inclusive amplitude' (including the 
bremsstrahlungsgluon) as an effective matrix element for b — > 57 as in equation (3.5) of 
||. There, this was just a mathematical trick to obtain a convenient form of the result 
up to the desired accuracy. Physically, we might envisage it as a replacement of the eight 
gluons by an abelian vector. This interpretation follows if one calculates the cross section 
of the bremsstrahlungsprocess which consists in summing over the final states and where 
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the sum over the gluons indeed yields an 'abelian' contribution (with the correct factor). 
In this way the final result is simply 



nindr^i C B CyJm B a s {m b ) row/, / / \ . a 
F incl [0 7 ] = - ^-f ^^(t?7 log(m 6 //x) + r 7 ) x 

dy-n -<M2/)0b - 

y/(l+y) V 2 / 

jpinclrn l C B Cy^/m B a s (mb) n . , , N 

F mct 8 = - — | J — ^ log m b //i + r 8 x 

l+y) V 2 / 



F-\0 2 ]= - CBC lf^ asi ™ b) (k log(mM + r 2 ) x 

^-77 -<t>v{y)<t>B 7, • 

yJ(l+V) V 2 / 

In (3) and (4) r§ and r2 include real and imaginary part ||. 

In Method 2 we cut the virtual loop momenta at a value A. The same cutoff A is then 
used for the momentum when calculating the Oj diagram (Fig. 3 a,b) with a gluon ex- 
changed between the quark lines (next section). No Bremsstrahlung corrections are then 
added. The physical picture here is that gluons with momenta below A are embodied in 
the hadronic wave function. The formulae for the contribution of the operator O7 with a 
cutoff are given in Appendix A. 



4 NLL "exclusive" amplitudes 

Next consider the 'exclusive' amplitudes, i.e. amplitudes which include the exchange of 
gluon between two quark lines. 

As mentioned above, we require the gluons to be of-shell by at least an amount A 2 . In 
the numerical evaluation, we will consider the three cases A =0.2,0.5,lGeV. 



4.1 Contributions from O2 

The form factors associated with O2 originate from two types of diagrams (Fig. 4a, b). 
When the photon is emitted from external quark lines (Fig. 4a) , we can use the expression 
for the effective bsg vertex given in ||, |14[ 



167r 2 ^ ' 2 \m 



4 



V 



rrii 



4r(e)/i e exp(7£;e)(l — e) exp(i7re) / dx[x{l — x)\ 



ll-e 



rri; 



x(l — x) 



+ iS 



The explicit calculation of V(r) (r 



gives: 



V(r) 
for r < 4, and 



2 1 2 
37 + 3 



-- + -log^ + -- + -- 

nij 9 r 3 r Jr{A-r) 



4 6 + r 4 (r 2 - 2r 



arctan 



4 - r 



(17) 



2 1 2 /i 2 4 6 + r 2 (r 2 - 2r - 8) - 

3 e 3 ° S m 2 9 r 3 r </ r (4 - r ) ° S + y/r — A 



2 r + 2 /r-4 

3 r V r 

for r > 4. We see that V(r) contains an infinity -, where e is the parameter which 
comes from the dimensional regularization method used. It is of course cancelled by an 
appropriate counter-term. The contribution from V(r) to the effective Hamiltonian then 
has the following form 



C 2 V{r] 



Ots_}_ 

An q 2 



(19) 



while that of the counter-term reads: 



3 48vre 

-~C 2 ^- 
3 48vre 



4 OL 

< su'-f\0 3 \bu > +3-C 2 — — < sujlOAbu > 
3 487re 
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< su-f\0 5 \bu > +3-CV 

3 487re 



< su^\Oo\bu > . 



(20) 



It is easy to see the cancellation of the divergences proportional to - using the formula 

2 Kib 2 ^cd = Q^ad^bc ~ &ab&cd)- 

Summing over all diagrams where the photon is emitted from different external quark 
lines we finally obtain for the form-factor in the rest frame of the B-meson 

9s N 2 -l 1 

- • ' " ■ " : 2 aJL, q L 



F 1 [0 2 



16tt 2 2N m\ 



(fi B (E q )(j) V (y) 



2(m B - E q ^ * 



El 

dE q dzdy x 



y i-yj 



(21) 



where Q u = 2/3, Q d = -1/3, n = ((1 - y)p v ~P q ) 2 , r 2 = ((1 —y)pv +P 7 ~Pq) 2 , z = cos9 
and 9 is the angle between the spectator quarks in the B and i^*-mesons, respectively. 
When the photon is emitted from the quarks in the internal loop (Fig 4b), we will use 
the following expression for the effective sbyj-vertex [11]: 

g s eQ u A a 



r 

fJ,U 



8tt 2 



iep^atq 13 Ai 5 +j£Ai 6 )+ 



p 7 ■ q ' P-y ■ Q 



(22) 
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where q is the momentum of the gluon. The quantities Ai n = Ai n (zo, -21,-22) are functions 
of the variables zq = s/m 2 ,zi = q 2 /m 2 ,Z2 = p 2 /m 2 , and s is the invariant mass of the 
internal quark pair. When the photon is on-shell we have 



Ai 5 (20,21,0) = -1 + 
Ai 6 (z ,z 1 ,0) = 1 + 



— (Qo(z ) - Qo(zi)) 

Zq — Z\ Zo — Z\ 

zi ,„ , , „ , „ 2 



(<2o(^o) - Qq{z\)) 



(Q-(zo)-Q-(*i)) 
(Q_(*b)-Q-(*i)) 



Zq — Z\ Zq — Z\ 

Ai 2 3(^o, zx, 0) = Ai 5 (z Q , z x , 0) = -Ai 2e (z Q , z x , 0). 



(23) 



The functions Qq and Q- are defined in the following way: 



-2iri log 



Qo(x) 



y/x + \Jx - 4" 



> (au-i)i«g« 



'£ -4 /- 1 , 

x du- 



The expression for the form factor is the following: 



f 2 [o 2 ] = a s c B c v c N J <p B {E q )My) 



1 — xu[l — u) 

u(2u- 1) 
1 — xu(l — u) 



E a 

dE q dzdy x 



2(m B - E q ) 16g 2 vr 3 



(24) 
(25) 

(26) 



x 



'4(1 - z)E q 
m B 



A i5 - 4A J6 + 2 



z 2 )E 2 



Pi ■ q 



(m b (l - y) - 2E q )(l + z)m b 
p 7 • q 



A 



'20 



4.2 Contributions from 0$ and O7 

The contribution of the diagrams where the photon is emitted from the spectator quark 
lines (Fig. 5) in the B and K mesons equals 



CbCvCnQu 



ft. 



47T 3 171% 



dE q dzdy - 



E q {m B - 2E q 



E q (2-y(l-z)) 



[l-y)^2{m B -E q ) 
-6 B (E q )(j) v (y). 



m B (l-y)-E q (2-y)(l-z 
Similarly, when the photon radiates from b- and s-quark lines (Fig. 5), one obtains 



-F 2 [Os] — C b CvCnQ(1 



47T 3 m 2 B 



dE q dzdy 



Eg{m B - 2E q 
yj2(m B - E q 



--(f) B {E q )(f) v {y). 



(27) 



:28) 



The contribution associated with the operator O7 (Fig. 3 a,b) has been calculated in 



H 



F[0 7 



16tt 4 



C N C B C V \ dydE q dz(p v (y)4> B (E q 



E n 



E 2 (m B -2E q )(2-y(l-z)) 



<r= {^-y) m K* -(m B {l- z)+m K *{l + z] 

m B v 

N b = 2E q m B - m|(l - y) - m 2 K *y(l - y). 



N b qy2(m B - Eq) 

(i-y) 



(29) 
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In general, we neglect the mass of the K* meson and of its constituents. However, 
when we deal with infrared divergent integrals and use the cutoff parameter A, we cannot 
neglect mj?. if it is larger than A. Therefore we take into account the finite values of rriK* 
and m s when calculating the matrix elements of O7 (Figs. 2 and 3). 

4.3 Exchange diagrams 

Exchange diagrams, which can contribute to the decays B —>■ V'j are represented in Fig. 
6 for the operators 0± u , C>2u, O3, O4, O5 and Oq. 

For B — > K* r y decay the contributions of the operators 0± u , 2u are proportional to the 
small CKM factor V u bV* s and consequently can be neglected. Also the contributions of the 
operators O3, O4 can be neglected in the approximation where the vector meson mass is 
equal to zero. The contribution from the operator O5 comes from the exchange diagrams 
shown in Fig. 6. The expression for corresponding form factor is the following: 

F[05l = -JL gggK J E M)Mv) (1 + ,) . ,30) 

The form factor for Oq is given by the same formula except for an additional factor 3 
coming from the trace in color space: i^[Oe] = 3-F[Os]. 



4.4 Renormalization scale dependence 

It is known that NLL contributions drastically reduce the large \x dependence of inclusive 
decay rate for B — > X s ^\ more exactly, the fi dependence of the LL term proportional to 
Cj is canceled by explicit logarithms, proportional to a s . As we have here an additional 
ji dependence in F 1 [0 2 ] and also in C 5 , C 6 in the exchange diagrams, it is necessary to 
check whether such a cancellation also materializes in the exclusive decays B — > Vj. The 
\i dependence of the contributions of the operators 5 , 6 is due to the \x dependence of 
the Wilson coefficients C5, Cq which is determined by the renormalization group equation: 

An approximative solution for C5, Cq is 

C, = ^7 2 ° 5 log (»/M w ) C 6 = ^7 2 ° 6 tog(li/M w ), 

where 72 5 =-2/9, 72 6 =2/3. The exchange diagrams contribute in the combination (36*6 + 
Ck)F[0§\. On the other hand the scale dependence originating from -F[0 2 ] is the same 
and opposite as seen from Eq. (|17D and the color matrix identity given there. 
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5 Numerical results and discussion 



The decay amplitude is given by the following expression: 

AG 

M(B -»• K*i) = —-=-em B e^v t (ie^ a pp"e v p^ + (p^v)pv^ ~ (p^Pv)e Vp ) 

x (c 7 F° [O r ] + C 7 F ind [Or] + C 2 F incl [0 2 ] + C 8 F incl [0 8 ] (32) 
+ C 2 F[0 2 ] + C 8 F[0 8 ] + C 7 F[0 7 ] + (C s + 3C 6 )F[0 5 ]) , 

with F[0 2 ] = F l {0 2 ] + F 2 [0 2 ], F[O s ] = F^Os] + F 2 [O s ]. 

We begin with a discussion of the individual contributions for one special case: method 
2 and A = 0.2, \i = 5GeV, fs = 160MeV p^j . The contribution associated with the oper- 
ator 7 consists of two pieces: an "inclusive" part (Fig. 2) and an "exclusive" part (Fig. 
3). The contribution of the first diagram (to the total rate) is -37%, while that of the 
second is -20% Q. Similarly, the contribution of the operator 2 is altogether 55%. 51% 
come from the "inclusive" diagrams (Figs. 3,4 of ||), while the exclusive contributions 
amount to 5% (Fig. 4a) and -1 (Fig. 4b). The operator 8 contributes only 2.3%, which 
consist of 3.8% for the "inclusive" part (Fig. 9 of |§) and of -1.5% for the "exclusive" 
part (Fig. 5). The contribution of the diagrams associated with the operators O5 and Oq 
is 6.2%. The 0(a a ) corrections to C 7 amount to -5%. The contributions of the additional 
factors in (10) and (12) equal -7%. 



Our numerical results for the B — > K*^ decay rate are given in Table 1 in LL and NLL 
approximation (for the two methods to regularize the IR divergences). The dependence 
on the cutoff A is shown in Fig. 7, while the variation with /# is given in Fig. 8. For 
fs = 160MeV and A=0.5GeV our prediction for the branching ratio is (5.18 ± 0.5) ■ 10~ 5 
(method 1) and (5.16 ± 0.5) ■ 10~ 5 (method 2) (the errors refer to the scale uncertainty). 
For both methods, they are in agreement with the experimental value BR(B — > K*j) = 
(4.2 ±0.8 ±0.6) x 10 -5 0, ^3[, although it is somewhat on the high side. We note that with 
the simple expression for the vector meson wave function 4>y{y) = 6y(l — y), the prediction 
for the branching is lowered by ~ 15%. The NLL contributions clearly reduce the scale 
dependence of the LL result. For the decay rate in LL approximation, the uncertainty 
due to the scale dependence is 27%. On the other hand, the NLL results vary between 4% 
and 11% for method 1 and between 4% and 21% for method 2 for lGeV > A > 0.2GW. 
This remaining scale dependence is larger than the one of the inclusive decay. However, 
the standard model yields indeed an 'unnaturally' low scale dependence for b — > S7 [25|. 
Therefore the variation with the scale seems reasonable to us. 

The uncertainty due to the unknown value of A is between 6% and 14% for method 1 
and between 10% and 26% for method 2. The results for the two methods for fi = 5GeV 
and 0.2GeV < A < lGeV differ no more than 14%. Thus in method 1, the theoretical 
uncertainty is quite small, in the vicinity of 20%. If this is a correct assessment, then it 
is clearly desirable to improve the experimental accuracy. 



1 Since we only calculate to order a s , all corrections we calculate contribute linearly. 
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The rather strong A dependence of method 2 is clearly a shortcoming of this work. 
It arises mainly from the logarithmic cutoff singularity of the gluonic corrections to the 
operator O7. A clear improvement would be to take into account three body Fock states 
(with an extra gluon) whose matrix element could remove the A divergence. This requires 
relations between two-and three-body wave-functions. In simple cases, such relations exist 
P^| . However, a realistic calculation is beyond the state-of-the art. Another source of A 



dependence lies in the "inclusive" contributions associated with the operators O2 and 0§. 
To isolate it, one would need a two- loop calculation with a finite cutoff; this is beyond 
the scope of this paper. We note here only that the corrections to F[02], F[Os] vanish as 
A goes to zero. We hope to return to these issues elsewhere. But if the three-body wave 
function indeed removes the low A singularity, then a relatively low value of A, say 0.5 
GeV, would reproduce best the correct value of the branching ratio. 

In conclusion, we have calculated the next-to-leading logarithmic corrections to the 
branching ratio of the exclusive B — > K*^ decay, using an "inclusive" method and an IR 
cutoff parameter A to remove possible infrared divergences. For fg = 160MeV and A=0.5 
GeV, the prediction for the branching ratio is (5.18 ± 0.5) • 10 -5 and (5.16 ± 0.5) • 10~ 5 
for two IR regularization methods we use. This corresponds to a 10% increase over the 
leading order calculation. The A dependence is sizeable, around 20%. We discuss possible 
ways to improve this uncertainty. 
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APPENDIX A 

In this Appendix we present the results for the matrix element of inclusive b — > 57 decay 
associated with the operator O7 when we use a cutoff (method 2). The corresponding 
diagrams for the wave function renormalization factors z| and and the vertex op- 
erator are shown in Fig. 2. Following ||, we introduce the renormalization scale in the 
form /1 2 exp(7£:)/(47r). Then the MS subtraction corresponds to subtracting poles in e in 
the dimensional regularization scheme. 

First we present the results for nonzero gluon mass and without cutoff. For that case 
Z2 is given by: 

i^wo^e/VY / \ r jd f 1 j 4m 6 -2(p'-g) 



E« = -^i^r y exp( e7E ) £Vg jf dx 



(q 2 — 2xp ■ q — r 



2\2 



dp 



P =m b 
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where //o is the gluon mass, r 2 = xml + (1 — x)^ — xp" 2 . Then the integral over q is 
performed and we obtain: 

A b = x 2 m 2 b + (1 - x)fi 2 (A.2) 

The expression for Z 2 ^ can be obtained from (A.2) via the substitution mj -> m s . The 
contribution of the vertex diagram Fig 2b is given by 

2a s em b f 4 f 1 AekR(p b ■ p s ) - 2(ekRqp s + p b qekR) 
V = 3^W J dq Jo dXl I dX2 ( g 2_ 2(p . g) _ t 2 )3 ( A - 3 ) 

where p — (x 1 - x 2 )p b + (1 - Xi)p s , t 2 = (x l - x 2 )(m b - pi) + (1 - x^im 2 . - p 2 ) + x 2 ^\. 
After integrating over q, V can be expressed in the following form: 

V =< S1 \0 7 \b > tree V, (A.4) 

VI = 2(Xs f 1 dx 1 r dx 2 x 2 {m 2 + m 2 ) 

37r Jo Jo x 2 /Iq + (1 - x 2 )((xi - x 2 )ml + (1 - xi)m 2 )' 

Finally, the contribution of O7 for the case of nonzero gluon mass is given by 

< S1 \0 7 \b > M0 =< S1 \0 7 \b > tree (l + V, + ^ log ^ + (yzfzf - l)) . (A.5) 

Now we proceed to the case when we introduce a cutoff parameter A. We calculate the 
integrals for the region \q 2 \ < A 2 and use then the previous results to obtain the integrals 
over the region \q 2 \ > A 2 . We have the following expressions: 

Ej? = ~^(2nr (f] e eMei E ) I d d q C dx , ^ = ^ f ^ 

?37T 3V 7 ^4^; ^ m7 7 k 2| <A 2 y 7 ( g 2 _ 2a; (p' . g) _ r 2)2 



fW - 1 = , (A.6) 



(^2 )k 2 |<A 2 dp p > 



2a s em b r 4 r 1 AekR(p b ■ p s ) - 2(ekRqp s + p b qekR) 

V ^<" 2 = 3^8^ L<^ d q Jo dXl L dX2 (g 2 -2( P - 9 )-t 2 )3 ^ 

We have the following integrals: 



J 2A 



^4 1 

a q- 



37T 3 J\q 2 \<A 2 (q 2 — 2x(p ■ q) — r 2 ) 2 

(A.8) 



S 2A ; = ^ / „. .dq-j^- 



37r 3 V| g 2 |<A 2 p' 2 (g 2 — 2x(p' • g) — r 2 ) 



2 
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(1) 



d A q 



1 



3tt 3 V| <A 2 (g 2 -2(p-g) - t 2 ) 3 



(A.9) 



V, 



(2) 



A 3vr 3 V|<a^ ^ 2 (g 2 -2(p-g)-t 2 ) 3 
The integrals over the momentum q are performed using a Wick rotation and we obtain: 



s (6i) = _ 
2A 67nr 2 p' 2 



- A 2 + r 2 - J A 4 + r 4 + 2A 2 (r 2 + 2x 2 p' 2 )- 



-2x 2 p' 2 



,(62) 
J 2A 



2 '2~ 



arctanh 1 + 



1 f 1 



2arp 



+ arctanh 



-2x 2 p 2 - r 2 - A 2 



A 4 + r 4 + 2A 2 (r 2 + 2xV 2 



, . "tj ( A 4 + r 2 (r 2 - 2x 2 p 2 ) - (A 2 + r 2 - 2x 2 p 2 WA 4 + r 4 + 2A 2 (r 2 + 2xV 2 ) + 
12txx a [p 



+2r 2 A 2 ) + Ax 



arctanh 1 + 



2x p 



2x 2 p' 2 - r 2 - A 2 



arctanh 



A 4 + r 4 + 2A 2 (r 2 + 2xV 2 / 



V, 



(i) 



A 2 (2p 2 + t 2 ) - t 4 W(A 2 + t 2 ) 2 + 4p 2 A 2 + t 2 ((A 4 + t 4 + A 2 (V + 2t 2 )) 



A 3vrp 2 (p 2 + t 2 )((A 2 + t 2 ) 2 + 4p 2 A 2 ) 

(A 2 (p 2 + t 2 ) + t 4 ),/(A 2 + t 2 ) 2 + VA 2 - A 4 (p 2 + t 2 ) - t 2 A 2 (3p 2 + 2t 2 ) - t 6 



(A.10) 



(2) 



A 



37rp 4 (p 2 + t 2 )^(A 2 + t 2 ) 2 + VA 2 
where we take t 2 = £2/^0, p 2 = (1 — X2){(x\ — X2)m 2 + (1 — xi)m 2 ), r 2 = (1 — Then 



[Z\ )| g 2| <A 2 and V\ q 2\ <A 2 are equal to (p 



Vi 



i g 2| <A 2 =< s^\0 7 \b >tree a s dx x J dx 2 (m 2 b + m 2 ) (yj^ - (1 - a^V^) 
[zf ) k 2, <A 2 = l + a s ( l dx ( 4- ((4m 6 - 2p')E 2 6 A 1) + 2p'S 



(62) 
2A 



(A.ll) 



1 + a. /* f 2m h _ 2S g) + 2n»* E g> + 25^) 

jo y dp dp J 



Then the expression for (Z^ )\ q ^\>A? is given by the difference between Zrp and (Z2 )| 9 2| <A 2 
while Vj g 2| >A 2 is given by the difference of V and Vj g 2| <A 2. For both cases, we can set the 
gluon mass equal to zero, i.e we take the limit /i — > 0. 
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Figure Captions 



Figure 1: Leading order contribution, associated with operator O7. 
Figure 2: Order a s corrections for the b and s quark wave function renormalization 
and the operator O7. 

Figure 3: Order a s corrections to the matrix element of the operator O7 with gluon 
exchange between the quark lines. 

Figure 4: Order a s corrections to the matrix element of O2 with gluon exchange 
between the quark lines and photon emission from external (a) and internal (b) quark 
lines. The crosses indicate the possible places of photon emission. 

Figure 5: Order a s corrections for the operator 8 with gluon exchange. The possible 
places of photon emission are labelled with a cross. 

Figure 6: The contributions of the operators Oj, i=l,2,3,4,5,6 with quark exchange. 
The possible places of photon emission are labelled with a cross. 

Figure 7: The branching ratio Br(b — > K* / y/10~ 5 as a function of log[l/A(Gel / )] for 
the two methods discussed and f B =lQ0MeV. 

Figure 8: The branching ratio Br(b — > K*^y/10~ 5 as a function of f B for the two 
methods and A = 0.2GeV 
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Table 1. 

The branching ratio BR(B — > i^*7)/10~ 5 in LL and NLL approximation for /z=2.5GeV, 
5GeV, lOGev, A = 1.0,0.5,0.2GeV and / B =160MeV. 





// = 2.5GeV 


/i = 5G e y 


/i = lOGeV 


Br LL (B -> K* 7 )/10- 5 


5.91 


4.65 


3.70 


Br NLL (B -> K* 7 )/10- 5 , Method 1,A =lGeV 


5.35 


5.59 


5.48 


Br NLL (B -> K* 7 )/10- 5 , Method 1, A =0.5GeV 


4.70 


5.18 


5.21 


Br NLL (B -> AT* 7)/l(T 5 , Method 1, A =0.2GeV 


4.55 


5.09 


5.16 


Br NLL {B -> A^J/IO -5 , Method 2, A =lGeV 


6.21 


5.97 


5.74 


Br NLL (B -> fsr* 7 )/10- 5 , Method 2 A =0.5GeV 


4.66 


5.16 


5.20 


Br NLL (B -> iT^/lO" 5 , Method 2, A =0.2GeV 


3.46 


4.40 


4.69 
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'0 7 a 



Figure 1: Leading order contribution, associated with operator O7. 
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Figure 2: Order a s corrections for the b and s quark wave function renormalization and 
the operator 7 . 
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Figure 3: Order a s corrections to the matrix element of the operator O7 with gluon 
exchange between the quark lines. 



0, 



0, 









u 1 


10 u 
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Figure 4: Order a s corrections to the matrix element of O2 with gluon exchange between 
the quark lines and photon emission from external (a) and internal (b) quark lines. The 
crosses indicate the possible places of photon emission. 
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Figure 5: Order a s corrections for the operator 0$ with gluon exchange. The possible 
places of photon emission are labelled with a cross. 



b s 




Figure 6: The contributions of the operators O;, i=l,2,3,4,5,6 with quark exchange. The 
possible places of photon emission are labelled with a cross. 
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Figure 7: The branching ratio Br(b — > K*^/10 5 as a function of log[l/A(Gel / )] for the 
two methods discussed and / B =160MeV. 
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Figure 8: The branching ratio Br(b — > K*^/10 5 as a function of f B for the two methods 
and A = 0.2GeV. 
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